Phonon Bloch oscillations in acoustic-cavity structures 
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We describe a semiconductor multilayer structure based in acoustic phonon cavities and achievable 
with MBE technology, designed to display acoustic phonon Bloch oscillations. We show that forward 
and backscattering Raman spectra give a direct measure of the created phononic Wannier-Stark 
ladder. We also discuss the use of femtosecond laser impulsions for the generation and direct 
probe of the induced phonon Bloch oscillations. We propose a gedanken experiment based in an 
integrated phonon source-structure-detector device, and we present calculations of pump and probe 
time dependent optical reflectivity that evidence temporal beatings in agreement with the Wannier- 
Stark ladder energy splitting. 
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Electronic Bloch oscillations, that is, oscillations of an 
electron induced by a constant electric field in the pre- 
cense of a periodic potential, 01 are a beautiful and clear 
example of quantum effects in solids. When an electric 
field is applied on a charged particle in a crystal, its 
wavevector increases with time. Thereafter, Bragg inter- 
ference leads to a velocity reduction, and finally to a sign 
change at the band edge. Notwithstanding its simplicity, 
for many years the issue was controversial and only quite 
recently the existence of electronic Bloch oscillations has 
been definitively established. In normal metals the large 
Brillouin zone and electron relaxation lead to an over- 
damped behavior characterized by Ohm's law. Instead, 
electronic Bloch oscillations are observable in SL's due to 
the Brillouin zone reduction. Q 

Very recently photon Bloch oscillation devices based in 
optical microcavities have been proposed, and related 
reflectivity and time-resolved optical transmission fea- 
tures have been observed in structures grown on porous 
silicon. The concept behind these devices is also quite 
nice and simple. The optical microcavities provide the 
discrete photonic states which can couple through pho- 
ton mirrors thus leading to photonic minibands. The 
spatial dependent energy gradient, equivalent to the elec- 
tric field for electrons, can be achieved, e.g., by varying 
the refractive index, [jj] As compared to electrons, photon 
dephasing mechanisms are less effective and thus Bloch 
oscillations are, in principle, more easy to observe. On 
the other hand, photonic minibands require a large num- 
ber of optical microcavities (around 30-50), amounting 
to total thicknesses larger than 40-50 microns, [s] This 
is too large for current molecular beam epitaxy (MBE) 
semiconductor technology. For this reason the reported 
structures were grown using electrochemical methods of 
porous silicon nanostructuring that do not have this limi- 
tation. 3 The drawback is that the quality of the samples 
is not as good, layer interfaces and refractive indices are 
not that well controlled, and optical residual absorption 
becomes an important issue in the performance of the 
structures. 

Sound in solid media is also described by a wave equa- 



tion, the relevant parameters being the material density 
and sound velocity, and the boundary conditions estab- 
lishing the connection between displacement and strain 
between two different materials. Extending the above 
concepts applied to electrons and photons to phononic 
structures, in this letter we propose semiconductor mul- 
tilayer structures capable of displaying acoustic phonon 
Bloch oscillations. The building block of the structures 
are acoustic phonon cavities as recently introduced by 
Trigo and coworkers. . The phonon wavelength is only 
a few nanometers (phonon frequencies being in the THz 
range), and thus the sample full size and required inter- 
face layer quality are achievable with actual MBE tech- 
nology. In addition, acoustic phonon mean free paths 
are large compared to the structure size (typically 
below a micron) and thus dephasing is not a critical is- 
sue. We present calculations of phonon reflectivity, dis- 
placement distribution, and time evolution upon exci- 
tation with a localized phonon source to illustrate the 
device behavior and to demonstrate the existence of a 
phonon Wannier-Stark ladder (WSL) and Bloch oscil- 
lations (BO). In addition, using a photoelastic model 
we calculate the Raman spectra which we show give a 
direct measure of the created phononic WSL. Phonon 
BO can be independently probed using coherent phonon 
generation techniques. For this purpose we propose a 
speciflcally designed device and we present calculations 
of time dependent reflectivity induced by femtosecond 
impulsions in this structure. 

A periodic stack of two materials with contrasting 
acoustic impedances reflects sound. The first fc = 
folded phonon mini-gap in a SL is maximum with the 
layers thickness (d) ratio given by di/vi = (i2/3w2, the 
stop-band and reflectivity of such a phonon mirror be- 
ing determined by the acoustic impedance mismatch 
Z — p\V\l P2V2 and the number of SL periods. 0, A 
phonon cavity can be constructed by enclosing between 
two SL's a spacer of thickness dc = m\c/2, where Ac 
is the acoustical phonon wavelength at the center of the 
phonon minigap. The cavity conflned modes corre- 
spond to discrete energy states within the phonon stop- 
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band, their width (i.e., the cavity Q-factor) being deter- 
mined by the phonon mirror reflectivity. |^ When a large 
series of phonon microcavities are coupled one after the 
other, the discrete confined energy states form phonon 
bands that resemble the minibands in electronic SL's. In 
order to display Bloch oscillations, the energy of the i-th 
cavity must differ from that of the (i — l)-th in a con- 
stant value. Such a linear dependence with position of 
the phonon cavity-mode energy, analogous to an electric 
field for electrons, can be obtained by tuning the cav- 
ity widths. Q We thus consider a multilayer structure 
where each unite cell consists of an acoustic phonon mir- 
ror made by (n-|- 1/2) periods of two materials with 
contrasting acoustic impedances (GaAs/AlAs in the ex- 
amples discussed here), followed by a A cavity (GaAs). 
This unit cell is repeated Nc times, with layer thicknesses 
increasing from the surface to the substrate so as to have 
a linear decrease of cavity-mode energy by steps of A. 

The stationary solutions of the acoustic waves in the 
proposed structure can be derived using a matrix method 
implementation of the elastic continuum model. |^ The 
calculations give («) the phonon field distribution in the 
different layers, (ii) the phonon reflectivity and/or trans- 
mission, and (iii) the variation along the structure of the 
energy bands associated to an infinite SL with the local 
unit cell. The latter are given by the condition — 1 < 
(flu -I- a22)/2 < 1, where an and 022 are the diagonal 
elements of the transfer matrix across each period of the 
structure. In Fig. Q we present results for a Nc — 25 
period acoustic cavity structure. The first unit cell is 
made by a 2.5 GaAs/AlAs period 59.3A/23.5A phonon 
mirror {n = 2) and a GaAs spacer tuned to 20 cm^^ 
(79A). The energy steps are given by A = 0.15 cm^^, 
and the structure is limited from the two sides by GaAs. 
Panels (a) to (c) display, respectively, the phonon band 
structure (black regions represent "forbidden" energies), 
the phonon reflectivity, and the phonon displacement dis- 
tribution as a function of position and energy. In the lat- 
ter panel, calculated for phonons entering from the left, 
darker regions indicate larger acoustic phonon intensities. 
Several features should be highlighted in these figures, 
(j) An acoustic phonon allowed band originated in the 
coupled discrete cavity modes is observed between two 
forbidden minigap bands (panel (a)). The energy of the 
bands decreases linearly with cavity number according to 
design. Three different spectral regions can be identified. 
Between 17.4 and 18.6 cm~^ (shaded in Fig.^, phonons 
are confined in a spatial region limited by the top and 
bottom of the lower and upper minigap bands, respec- 
tively. On the other hand, below(above) 17.4(18.6) cm~^ 
a phonon entering from cavity #1 will bounce back at the 
bottom of the lower (upper) minigap band and leave the 
sample, (ii) The phonon reflectivity (panel (b)) displays 
a broad stop-band (between ~ 15 and ^ 21.5 cm~^), ba- 
sically determined by the superposition of the individual 
cavity minigaps. A series of dips modulate the reflectiv- 



ity within this wide stop-band. Between the lower stop- 
band limit and ~ 17.4 cm~^, and between ~ 18.6 cm~^ 
and the upper stop-band limit, a series of features with 
varying energy spacing can be identified. These originate 
from phonon interferences determined by a propagation 
limited by the sample surface and by the minigap bands. 
These features are relatively weak and broad because of 
the small acoustic impedance mismatch between the top 
GaAs layer and the cavity structure. On the other hand, 
sharp reflection dips are observed in the region between 
^ 17 A and ~ 18.6 cm~^. These dips, which are equidis- 
tant and separated by A = 0.15 cm^^ correspond to 
the coupling of external phonons, by tunneling through 
the minigap band, to the Wannier-Stark states confined 
within the structure. (Hi) This phononic WSL can be 
also identified by the well defined discrete phonon modes 
displayed in panel (c). The WSL is the spectral domain 
counterpart of the BO. It is precisely in this spectral re- 
gion where oscillations should appear in the time-domain. 

The time and spatial variation of the acoustic phonon 
displacement field Ug{z,t), created by a pulse described 
by a spectral function g{uj) and incident at t — 
at the GaAs-sample interface {z = 0), can be evalu- 
ated using the scattering method described by Malpuech 
and Kavokin. 01 Within this description, Ug{z,t) = 
^ u{z)g(uj)exp{—iujt)doj, where u{z) are the sta- 
tionary solutions of the elastic wave equation with fre- 
quency uj shown in Fig.^c). The time evolution of such 
wavepackets for the two different energies indicated in 
Fig. n^a) are shown in Fig. [3 For energies above or be- 
low the phonon WSL region (20 cm in the example 
shown), the incident pulse propagates within the sample 
up to a position where it is back-reflected by a mini- 
gap band leaving afterwards the sample. On the con- 
trary, when UJ corresponds to the WSL energy region 
(18.3 cm in the displayed figure), a fraction of the 
pulse energy is backreflected at the surface while another 
part enters the structure by tunneling through the mini- 
gap band developing afterwards clear periodic oscillations 
within the structure. In order for these phonon BO to be 
observed, the FWHM of g{uj) should be larger than A. 
For these calculations we have used a gaussian distribu- 
tion with 2(7 = 1.0 cm""'^. On the other hand, the period 
of the oscillations (and consequently the length travelled 
by the pulse) is directly determined by A {tb — h/ A). 

In what follows we will show how Raman scattering 
and coherent phonon generation experiments can pro- 
vide a direct probe of the phonon WSL and BO, respec- 
tively, in these devices. Raman scattering has extensively 
shown to be a powerful tool to study phonons in semicon- 
ductor multilayers and, in particular, confined modes 
in acoustic cavities. 0, In order to evaluate the Ra- 
man spectra, we use a photoelastic model. 0, 0| We an- 
alyze two experimental geometries, namely backscatter- 
ing (BS, ks ~ —kL, q =^ ^k^) and forward scatter- 
ing (FS, ks ~ k^, q ^ 0). Here k refers either to the 
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scattered or laser wavevector, and q is the transferred 
wavenumber. In Fig. 3 we present BS and FS spectra 
calculated for the sample described above and assuming 
laser excitation with 550nm. For comparison purposes 
the corresponding phonon reflectivity is also shown. The 
Raman spectra display a complex series of peaks in the 
stop-band spectral region. We have verified that such 
rich spectra are a kind of sample finger-print that can be 
used as a characterization tool. Interestingly, the BS and 
FS spectra display clear peaks and dips, respectively, at 
exactly the WSL energies. High resolution Raman set- 
ups working in the visible can discern spectral features 
with resolution better than 0.02 cm~^, [li| thus provid- 
ing a spectral-domain tool able to probe the underlying 
phonon WSL. 

Coherent phonon generation is termed the impulsive 
generation of phonons using high power ultrashort laser 
pulses. In the case of THz vibrations, femtosecond 
pulses are required. To the best of our knowledge the gen- 
eration mechanism is still an open issue, and no complete 
theory is available to describe the processes involved. We 
briefly describe next a model for pump and probe coher- 
ent phonon generation and detection based on a photoe- 
lastic coupling between light and acoustic phonons. 
This mechanism is the only active when pump and probe 
are below the gap. Above the gap other mechanisms can 
contribute, but we expect the main conclusions to re- 
main essentialy valid. Any arbitrary time and position 
dependent phonon displacement in the structure w{z,t) 
can be expressed in terms of the phonon eigenstates as 
w{z, t) = J r^u{z)sin{ujt)djjj. Assuming that phonons 
are generated coherently through a photoelastic mech- 
anism by a femtosecond pulse (modelled as Eo{z,t) oc 
d{t)exp{ikLz)), the coefficients for the above expansion 
can be obtained as — ^ J^^ p{z) \Eo{z)f ^"^^^^ dz. 
Here L is the length of the sample, p{z) is the photoelas- 
tic constant which is assumed constant in each layer, and 
du{z)/dz is the strain associated with an eigenstate of 
energy lo. We note that the FS {q = 0) Raman cross sec- 
tion discussed above is proportional to |r(j|^.|3| Once this 
excitation is generated, it evolves according to the time 
dependence sin{Ljt) and can be detected by a delayed, 
lower power, probe pulse that senses the time depen- 
dent change of reflectivity. The change in reflectivity 
can be calculated as Ar(t) oc Ae{z,t)exp{2ikz)dz, 
where the probe pulse has been assumed to be propor- 
tional to e''^^. [13 Introducing again the photoelastic 
coefficient to relate Ae with the strain dw{z,t)/dz, the 
time dependent change in reflectivity can be written as 
Ar{t) oc /q°° p{z) ^^qI'*^ exp{2ikz)dz. If ^(z) is constant, 
this equation implies that only phonons with q ~ 2k can 
be probed. On the other hand, the detector's p{z) can 
be designed to access other excitations by backfolding 
the phonon dispersion. Comparing this equation with 
the Raman cross section, |5j] it is easy to see that the 



observable phonon spectrum is related to the BS Raman 
scattering component of the detector structure. 

In order to generate the quasi-monoenergetic phonons 
required to induce BO, and to monitor the time evolu- 
tion of the latter, we have conceived a monolithic source- 
sample-detector device (see the scheme in Fig. 0J. The 
first GaAs/AlAs SL (SL/s in Fig. ^ acts as the phonon 
source and is designed to generate, by excitation with a fs 
impulsion above the gap, an elastic pulse with energy cen- 
tered at 18.3 cm~^ and width equal to ~ 0.8 cm~^. |ll| It 
is made of 30 periods of 43.21/51.51 GaAs/AlAs. The 
layer widths determine the energy of the generated pulse, 
while the number of periods define, due to finite size ef- 
fects, its width. The coefficient r^j of the phonon pulse 
(or equivalently the FS Raman spectra) generated in this 
SL is shown in the inset on Fig. 01 Once generated, the 
coherent phonons propagate into the structure and act 
as the phonon source g(uj) for exciting Bloch oscillations. 
Between the SL source and the structure, a GaAs 200nm 
buffer layer is introduced to screen the residual pump 
power, not absorved in the SL, from impinging into the 
cavity structure and generating unwanted frequencies. 
The g(uj) phonon pulse, on the other hand, propagates 
through the GaAs layer basically unaltered. The cavity 
structure is identical to the one described above. Once 
within this structure, Bloch oscillations develop and part 
of the energy is lost to the substrate upon each turn. 
Their effect at the right end side of the sample is probed 
by the second SL (SL/d in Fig.^, which acts as an en- 
ergy selective detector of the Bloch oscillations through 
their effect in the time variation of the reflectivity. To 
keep this time variation simple, a second GaAs layer is 
introduced to stop the probe beam from being modulated 
also by the cavity structure. The BS Raman spectrum of 
the SL defines the detector's bandwidth, which we have 
designed to include the pulse g{uj)- It consists of a 20 
periods 38.ll/68.ll GaAs/AlAs SL. Its BS spectrum, 
calculated for a 750 nm probe pulse, is shown in the inset 
on Fig. 21 The device is terminated by a thick (> 1/xm) 
Gao.5Alo.5As layer which acts both as a stop layer for 
chemically etching the GaAs substrate, and as a window 
for accessing the detector SL from the back. In Fig. 01 we 
present calculations of the reflectivity change as a func- 
tion of time. Fast oscillations dominate the reflectivity, 
corresponding to the frequency of the coherently excited 
phonons (18.3 cm~^). In addition, this fast component 
is amplitude modulated by an envelope whose frequency 
is determined precisely by the Bloch oscillation period. 
The Fourier transform of the time dependent reflectivity 
clearly shows the WSL frequencies within an envelope 
determined by the input source gij-^j) and the detector 
bandwidth. 

In conclusion we have extended concepts previously 
discussed in the context of electronic and photonic prop- 
erties of solids to acoustic phonon physics. We have 
described semiconductor structures based in recently 
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reported acoustic cavities and achievable with actual 
growth technologies, displaying phononic Stark ladders 
and capable of sustaining Bloch oscillations. We have 
also shown how Raman scattering and coherent phonon 
generation provide spectral and time domain probes, re- 
spectively, of these acoustic phenomena. Structures as 
the one described here can be exploited to enhance the 
coupling between sound and other excitations (electrons 
and photons), and as the feedback high-Q resonator of 
a phonon laser. Moreover, engineered phonon potentials 
are not limited to linear dependencies that mimic electric 
fields but can take arbitrary shapes. This opens the way 
to novel phonon devices based in the discrete confined 
states of acoustic cavities. 
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FIG. 1: (a) Phonon band structure. The central white region 
hmited by the thicker black forbidden bands corresponds to 
the cavity-mode miniband. The white and thinner black re- 
gions above and below this bands arise from Bragg oscilla- 
tions. The arrows indicate the energy of the pulses displayed 
in Fig. 2. (b) Phonon reflectivity. In (a) and (b) shaded areas 
indicate the WSL region, (c) Phonon displacement as a func- 
tion of position and energy. Darker regions indicate larger 
acoustic phonon intensities. The arrows label the WSL ener- 
gies. The solid lines limit the "forbidden" bands. See text for 
details. 



FIG. 2: log{u{z, t)) for gaussian phonon wavepackets centered 
at the two energies indicated with arrows in Fig. H^a) , and 
FWHM given by 2cr = 1 cm"\ 



FIG. 3: Back and forward 550 nm Raman scattering spectra 
corresponding to the acoustic cavity structure with Nc = 25 
described in the text. For comparison the corresponding 
phonon reflectivity is also shown. The inset highlights the 
Stark ladder spectral region (17-19.5 cm~^). 



FIG. 4: Top-left: Source and detector SL's Raman spectra. 
Top-right: Scheme of the proposed monolithic source-sample- 
detector phonon device (see text for details). Bottom-left: 
Time variation of the reflectivity probed at the detector SL. 
Bottom-right: Fourier transform of the detected reflectivity 
change. 
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